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Production mechanisms for quarkonium states in hadronic collisions remain difficult to under- 
stand. The decay angular distributions of J/ip or T(nS) states into fi + \i~ final states are sensitive 
to the matrix elements in the production process and provide a unique tool to evaluate different 
models. This talk will focus on new results for the spin alignment of T(nS) states produced in pp 
collisions at ^/s = f .96 TeV using the CDF II detector at the Fermilab Tevatron. The data sam- 
ple corresponds to an integrated luminosity of 6.7 fb -1 . The angular distributions are analyzed as 
functions of the transverse momentum of the dimuon final state in both the Collins-Soper and the s- 
channel helicity frames using a unique data-driven background determination method. Consistency 
of the analysis is checked by comparing frame-invariant quantities derived from parametrizations of 
the angular distributions measured in each choice of reference frame. This analysis is the first to 
quantify the complete three-dimensional angular distribution of T(IS'), Y(2S) and T(3S) decays. 
The decays are nearly isotropic in all frames, even when produced with large transverse momentum. 



I. INTRODUCTION 

The physical process by which a heavy quark- 
antiquark pair is created in a hadronic collision and 
then emerges bound as a colorless hadronic state 
has been the subject of intense discussion within 
the framework of QCD for many years. The change 
of color configurations by soft gluon emission was 
first described in the Color Singlet Model (CSM) by 
Baier and Ruckl [l| in 1983. Subsequent cross section 
measurements at CDF I |2| in 1994 showed that the 
predictions underestimated the data by factors of 
10-50. An innovative effective field theory approach 
(NRQCD) by Bodwin, Braaten, and Lepage [3] recast 
the problem in terms of a perturbative treatment at 
the quark level with long-range process-independent 
parameters to describe the hadronization. The 
free parameters in the problem were set from the 
differential cross section data but result in a definite 
statement about the contributing matrix elements in 
the process. Therefore, the polarization state of the 
resulting vector meson is determined from the theory. 
The prediction is that at large px the quarkonium 
state should be transversely polarized. 

The prediction refers to quarkonium states that 
are directly produced in the hard collision. Experi- 
mentally, one measures quarkonium states that are 
prompt, i.e., that come from the primary production 
vertex. Therefore, the experimental data include not 
only directly-produced quarkonium states but also 
those that come from the decay of higher-lying QQ 
excitations, either S-wave or P-wave, a process known 
as feed-down decays. These extra events will almost 
certainly have different polarizations than the directly 
produced states. This presents a fundamental diffi- 
culty in interpreting the results from all experiments 
thus far. Some theoretical calculations attempt to 



estimate the feed-down depolarization effects, but 
there are large uncertainties in these estimates. 

Recent phenomenological approaches have taken 
different viewpoints than the NRQCD approach 
to describe the quarkonium production process. 
In the kx factorization approach of Baranov and 
Zotov 4], the color-singlet process dominates, and 
the heavy quark pair is isolated from the rest of the 
interaction by a minimum momentum transfer from 
each parent hadron. They predict a rather featureless 
Pt dependence that should be longitudinal for direct 
production, as will be seen later. In contrast, the 
higher-order CSM approach of Artoisenet, et al. @ 
considers multiple-gluon interactions and finds that 
the higher-order processes contribute coherently to 
the production amplitude and so build up the cross 
section to a point that generally agrees with current 
quarkonium observations. This is a dramatic shift 
and has reinvigorated interest in this model. 

In all three cases there are unique predictions 
for the polarization expected in direct production 
of quarkonium. This makes these new Tevatron 
results on T(nS) polarization and the expected new 
LHC results on quarkonium polarization especially 
important. 



II. PREVIOUS QUARKONIUM 
POLARIZATION STUDIES 



A. Charmonium Polarization Results 

Following the success of NRQCD predictions 
for the J/ip and ^(2S) differential cross section 
results, using the freedom available in matching the 
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long- distance parameters of the model, CDF set 
out to improve the measurements of the prompt 
polarization for these charmonium states using the 
larger event samples of Tevatron Run II and the 
improved measurement accuracy of the CDF II 
detector. The prompt polarization for both J /if} and 
-0(2S) does not increase as pr becomes significantly 
larger than the mass of the state, i.e., in the pertur- 
bative regime. The results for the J/ip are shown 
in Fig. [TJ The measured polarization is small at all 
Pt values, in contrast with NRQCD expectations. 
The for prediction for direct polarization has a 
large magnitude, but the shape matches experiment. 
Depolarizing effects from feed-down might reduce the 
discrepancy in size. More work is needed in this area. 
These results are the state of the art at this moment. 
New results are expected from LHC experiments soon. 
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FIG. 1: Prompt J/ip Polarization Parameter a in the s- 
channel helicity frame as a function of the dimuon vt- 
Theory curves for NRQCD and kr factorization are in- 
cluded for comparison. The NRQCD band includes esti- 
mates of feed-down depolarization effects. 



B. Bottomonium Polarization Results 

The charm quark mass is close to Xqcd, and there 
is some speculation that disagreements with NRQCD 
predictions might imply that the charmonium system 
is too light for NRQCD assumptions to be valid. That 
is certainly not the case with the bottomonium sys- 
tem, which is the focus of the rest of this report. An- 
other issue whose importance has been emphasized in 
recent years is the choice of quantization axis along 
which to report polarization. In fixed target exper- 
iments the Collins-Soper (CS) axis [6(, the bisector 
of the angle between the colliding hadrons in the rest 
frame of the dimuon system, was the axis of choice. In 
hadron colliders a natural axis was the s-channel helic- 
ity axis (SH), defined as the direction of the Lorentz 
boost that transforms the dimuon from the lab sys- 
tem to its rest frame. In elementary quantum me- 
chanics we learn that a fully-polarized electron will 



have no spin alignment along an axis perpendicular 
to its quantization axis. Braaten et al. [7[ emphasized 
the need to analyze experiments in several different 
Lorentz frames to ensure that a small polarization 
was not an accident of frame choice. More recently 
Faccioli et al. Q have emphasized the importance of 
measuring the three independent parameters that de- 
scribe the parity-conserving decay of a vector meson 
into a pair of spin- 1/2 fermions in the vector meson 
rest frame: 

cx — [1 + \gcos 2 8 + X^sin^O cos2<j) 

ail 3 + Ag 

+ Xg c j > sin20 coscf)] 

The new results reported here will follow the 
admonitions from these papers and constitute the 
first hadron collider results that measure all three 
polarization parameters in each of two Lorentz 
frames: CS and SH. In addition Ref. [8[ pointed out 
a frame-independent quantity that serves both as a 
systematic constraint on the measurements in two 
or more Lorentz frames and also, should it be small, 
indicates that there is no significant polarization 
in any Lorentz frame. We shall also discuss that 
quantity for the new CDF results. 

Just as in the charmonium system, all three pro- 
duction models (NRQCD, for, and higher-order CSM) 
contain enough free parameters to match the T(15) 
cross section data well over the whole measured px 
range. The polarization measurement once again 
stands as the essential test of validity for the model 
results. Prior to this new CDF Run II measurement, 
there were results from CDF Run I [9( and DO Run 
II both in the SH frame. The two experiments 
do not agree with each other, and neither experiment 
matches theory, as can be seen in Fig. [5] The ex- 
perimental disagreement is not unique in polarization 
measurements. These are differential measurements, 
and good control of acceptance and efficiency is cru- 
cial to obtaining stable results. This experimental fact 
of life makes using the frame-independent test an im- 
portant validation of new results. 

One major difficulty with the T(nS) measurements 
is handling the dimuon background, which is much 
larger than in the J/ip case. The CDF and DO back- 
ground problems are very different, as can be seen in 
Figs. Eland Fig. H 

While it is clear from the preceding two pictures 
that the background fraction under any T(nS) peak 
is smaller for CDF than for DO, the importance of 
that fact really comes from the CDF observation that 
the combinatoric background has an angular variation 
that changes dramatically with dimuon mass and px- 
At the Tevatron the dominant dimuon background 
comes from heavy flavor decays. CDF defines a dis- 
placed dimuon sample by selecting dimuon-triggcrcd 
events in which at least one muon has a miss-distance 
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FIG. 2: T(15) polarization parameters as a function of pt 
prior to the new CDF measurement. The disagreement 
between CDF Run I and DO Run II is dramatic and is 
still not understood. The two theoretical curves for kr 
represent extreme choices for the feed-down effect. The 
CSM curve does not include feed-down. 



with respect to the primary vertex that exceeds 150 
/Mii. This selection defines a sample dominated by 
combinatoric dimuons, with a tiny admixture of mis- 
measured T(nS) events. We divide the dimuon mass 
region into twelve intervals as shown in Fig. [3J three 
regions containing T(nS) signals plus background and 
nine pure background regions. We make the funda- 
mental assumption that for the combinatoric dimuons 
the angular distribution is independent of the appar- 
ent decay length, i.e., the uncorrelated dimuons can 
appear to have a vertex at a distance that is un- 
correlated with their angular characteristics in their 
rest frame. If this hypothesis is true, then the angu- 
lar characteristics of dimuons in the pure background 
mass slices will be the same for the prompt and dis- 
placed samples. We have tested this hypothesis in all 
nine background bins and obtained excellent agree- 
ment. We show as examples the normalized distribu- 
tions in cos Ocs and <fics for two such bins, one at low 
mass, one at high mass, in Fig.[5j We make these com- 
parisons for all three parameters in both CS and SH 
frames and see comparably good matches. This agree- 
ment shows that we can extract the angular charac- 
ter of the combinatoric background under each of the 
three prompt T(nS) signal regions by using the mea- 
sured displaced sample in that same mass region. No 
sideband extrapolation is needed. 

To establish the normalization for the displaced 
sample, we plot the ratio of the two dimuon samples 
as a function of dimuon mass, shown in Fig. [5] The 
ratio is essentially constant in the nine background- 
only samples. We fit the ratio to a linear function 
of dimuon mass and use that fit to scale the displaced 
sample angular distribution to predict the background 
prompt angular distribution in the three signal bands. 

We fit the background angular distributions to the 
same functional form as the vector meson signal but 
do not limit the three coefficients for the background 



to the range expected for a real vector meson. We 
make a simultaneous fit to the displaced angular 
distribution and the prompt angular distribution 
in each of the twelve dimuon mass ranges. The 
prompt angular distribution contains a signal term (if 
appropriate) and an extrapolated background term, 
using the ratio as shown above. All fits give good x 2 
results in both CS and SH frames. We also included a 
cos A background term. This term cannot contribute 
for a real vector meson decay and it plays no role in 
the polarization analysis. 

An example of the polarization fits is shown in Fig. [7] 
for the lowest pt bin. The twelve dimuon background 
regions are each represented by a point in the two- 
dimensional parameter space in Xg, A^, along with the 
three signal fits. Error ellipses for la and 2a uncer- 
tainties are shown for each point. As can be seen 
from the plot, the background polarization parameter 
Xg changes non-uniformly with dimuon mass in this 
Pt region. Thus, having an accurate representation 
of background angular behavior from the data itself is 
essential to extracting the correct signal polarization. 
This capability is unique to this analysis. 

To establish systematic consistency for the measure- 
ment, one can use the frame- independent quantity 
A = -yzx~^ defined in Ref . [1] to compare results from 
the CS and SH frames. Because the same data are 
used in each frame, one has to evaluate matching like- 
lihoods by using a toy Monte Carlo procedure. In all 
Pt bins the results from the two frames agree well, 
based on the Monte Carlo expectations. 



III. RESULTS 

The complete results from this analysis, includ- 
ing tables of all fitted parameters in the CS and 
SH frames along with the frame-independent param- 
eters are published in the supplemental material to 
Ref. [il) . We give illustrative examples here, starting 
with those for T(1S) prompt decays, shown in Fig. [8] 

The T(1S) polarization parameters are all small 
and show little variation with pt in either the CS 
or SH frames. The values of the frame independent 
parameters are also small and p-r-independent, 
confirming that there is no Lorentz frame in which 
one would see large polarization for this state. 

To avoid having to interpret depolarizing effects 
of feed-down to T(15), we can look at the T(35) 
parameters. There is much less feed-down to this 
state because of smaller branching fractions from the 
X3b states and because the higher- lying T(nS) states 
have very tiny branching fractions to T(3S*). This is 
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FIG. 3: CDF dimuon mass spectrum integrated over pr with | y(fifi) |< 0.6. Note that the three T(nS) peaks are well 
separated. The different vertical bands denote dimuon mass slices used to study signal and background. 
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FIG. 4: DO dimuon mass spectrum for low and mid-range 
Pt bins integrated over the rapidity range | y{^^) |< 1.8. 
The three T(nS) peaks have pr-dependent overlap. 



the first measurement of polarization parameters for 
hadroproriuction of the T(3S) state. The sample size 
is smaller than that for the T(1S), but the results 
shown in Fig. |H] give no indication of significant po- 
larization at high pt in either the CS or SH frames. 

The new CDF results for the T(15 l ) agree well with 
the CDF Run I measurement and disagree strongly 
with the DO measurements, as shown in Fig. [TU1 
The situation for the polarization parameter Xg, 
also known as a, for the SH frame is summarized 
in Fig. 1101 There is no indication of a trend toward 
NRQCD behavior, illustrated by the yellow band. 
Shown on the plot are two curves for the kx model 
including feed-down from Ref The solid line 

assumes spin coherence between Y states, while the 
dashed line shows the case when the feed-down T(IS') 
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FIG. 5: Normalized distribution of the CS helicity polar- 
ization angle 9 and azimuth (j> for dimuons from the prompt 
dimuon sample (solid histogram) and the detached dimuon 
sample (points) for low mass and high mass regions that 
are pure background. The agreement in shape validates 
the claim that the angular characteristics of the combina- 
toric background are independent of the apparent vertex 
displacement. 



is depolarized. The data seem to favor depolarization. 
A direct measurement of T(1S) polarization from Xbi 
parents is clearly quite important for understanding 
the situation. Not shown here but included in Ref. [5( 
is a polarization band for the higher-order CSM that 
has the same general shape as the data and is consis- 



FPCP2012-12 



Presented at Flavor Physics and CP Violation (FPCP 2012), Hefei, China, May 21-25, 2012 



5 



~- < 


i 


I 


CMUPi 


CMU, 4-.0<p T <6.0G»Vfc 

i 3? 









order CSM model but not NRQCD expectations. We 
await the upcoming analyses from LHC to see the 
trends at higher energies and higher px in pp colli- 
sions. 



FIG. 6: Ratio of dimuon mass distribution for prompt 
and displaced dimuon samples. The ratio is essentially 
flat in the background regions. We use this ratio to scale 
the displaced vertex background angular distribution in 
the three T(nS) signal bands. No fit to the signal mass 
function is needed. 

tent with these results. The width of the uncertainty 
band can be reduced by having more experimental 
information about feed-down contributions. 



IV. CONCLUSIONS 

These new CDF results represent the first hadron 
collider determination of the complete T(15),T(25) 
and T(35) polarization parameter set in both Collins- 
Soper and s-channel helicity frames for promptly- 
produced bottomonium decays to /i + ' \i~ . The results 
show that all polarization parameters are small and 
have little variation with dimuon px- The trend of the 
data tends to support the kx model and the higher- 
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FIG. 7: Angular fit parameters Xg vs A^ for three signal and twelve background prompt dimuon mass regions for 2 
GeV/c < pr(i-ii-i) < 4 GeV/c. The signal fits cluster near zero for both parameters, while the background fits tend to 
have sizeable (f> variation and a mass-dependent range of 6 variation in the s-channel helicity frame shown here. 
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FIG. 8: Fit results for the three polarization parameters Ae,A^ and Xecf, for the Y(IS') signal as a function of pr- In 
neither the CS or SH frames is there any significant polarization effect. Systematic uncertainties are shown as extensions 
of the statistical error bars on each point. 
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FIG. 9: Fit results for the three polarization parameters Ae,A^, and Xe^, for the T(3S) signal as a function of pr- In 
neither the CS or SH frames is there any significant polarization effect. Systematic uncertainties are shown as extensions 
of the statistical error bars on each point. 
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FIG. 10: T(IS') polarization parameters as a function of pr for the new CDF measurement. The disagreement with DO 
Run II persists. The theoretical curves for the kr model is given for two possible feed-down situations, as discussed in 
Ref. 0| . 
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